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I. INTRODUCTION 


The experimental lower limit on the lifetime of the proton is a severe problem for models 
in which the fundamental scale of quantum gravity is low compared to the supersymmetric 
GUT scale 10^® GeV [^0. A baryon number U(l) symmetry cannot be gauged in field 
theory; like other accidental symmetries it is expected to be violated by effects at the string 
scale or by quantum wormholes and virtual black holes [Q]. Such B violation would appear 
at low energies as nonrenormalizable operators, for example of the form Xqqql/Mj where 
Mf is the fundamental scale and A is expected to be 0{1) in the absence of suppression 
mechanisms. Then for Mf in the 0(10 — 100) TeV range, for which collider signals of 
the fundamental degrees of freedom or of large extra dimensions may be observable, Tp 
comes out to be under a second, to be compared with the (mode-dependent) experimental 
bounds of order 10^^ years Various solutions have been proposed 1 , 1,0 all of which 
have implications for the production of an excess of baryons over antibaryons in the early 
Universe, for which B violation is a precondition. 

The production of gravitational Kaluza-Klein modes in the early Universe in such mod¬ 
els, in which some compactified extra dimensions are orders of magnitude larger than the 
fundamental length, gives a severe upper bound on the maximum temperature attained con¬ 
sistent with cosmological observations. Even for the maximum number (usually considered 
as 6) of large extra dimensions and the relatively large value Mf = 100 TeV, a maximum 
temperature of a few GeV is the upper bound if overclosure of the Universe, disruption of 
the successful predictions of nucleosynthesis, and an observationally unacceptable level of 
background gamma-rays from K-K mode decay are to be avoided |^. The energy density 
in K-K modes would also affect the evolution of density fluctuations by altering the time of 
matter-radiation equality Astrophysical production and decay of such modes also leads 


to an independent lower bound on the fundamental scale [10|, which is also constrained by 


the non-observation (to date) of direct and loop effects in experiments |n 


Any attempt at explaining proton longevity and baryogenesis should operate within these 
constraints. Exact (anomaly-free) discrete or horizontal symmetries can be imposed to forbid 
5-violating operators mediating proton decay |^, while allowing others, through which 
baryogenesis occurs: this approach requires an “X-boson” to be present, with couplings 
which appear unnaturally small (in contrast to the standard GUT or leptogenesis scenarios). 
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Baryon number can be gauged if the anomaly is canceled by a string theory mechanism, 
or B violation may be forbidden to all orders in perturbation theory by string selection 
rules, in some “intersecting brane” models |^. Note however that in a more general class of 
intersecting brane models 0, such selection rules do not prevent the four-fermion operators 
from appearing, as discussed above, in which case the fundamental scale cannot be low. 

If baryon number is perturbatively exact, nonperturbative processes [|H| are the only 
option to create net baryon number. It is difficult to see how this proposal can be recon¬ 
ciled with cosmological constraints, since any such processes would operate at or above the 
electroweak scale and be enormously suppressed at low temperature . 


A novel mechanism for suppressing proton decay in the context of a low fundamental 
scale theory has been proposed in Ref. 0 . This is a geometrical mechanism for suppressing 
4D R-violating operators, namely localization of fermions in extra dimensions. The simplest 
implementation is for the SU( 2 ) 2 , x U(l)y gauge helds to propagate in one extra dimension 
(c/. [0]), in which the quark and lepton wavefunctions are peaked about points separated 
by a distance L 30Mj ^ (Fig. !)• Then any strong R-violating operators in the effective 
5D theory can only produce proton decay proportional to the overlap of the wavefunctions, 
which can be exponentially small. Alternatively, proton decay by exchange of massive 
modes is suppressed by the Yukawa propagator over the distance L |^. Nonperturbative 
quantum effects which may lead to proton decay, for example virtual black holes []^, are 
also exponentially suppressed due to the integration over the hfth dimension. 


Baryogenesis in such scenarios has been discussed in Ref. |^. However, the authors of 
that work introduced a Standard Model singlet scalar with renormalizable baryon number- 
violating interactions and primarily considered the effect of thermal corrections to the wave 
function overlap which controls the strength of the baryon number violation. 

An Affieck-Dine-type scenario consistent with large extra dimensions, which involves 


introducing a scalar charged under a new U{1) group, has recently been described in |]T8 


while the idea appears viable from a cosmological viewpoint, the authors do not discuss how 
R-violating couplings which induce proton decay are to be suppressed, while still keeping the 
operators by which the scalar decays. One might imagine that additional exact symmetries 
could be imposed to achieve this, in the spirit of p. 

In this work, we consider a related, but simpler class of scenarios and report mostly on a 
“no go” result. Specihcally, we consider the feasibility of baryogenesis under the following 
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conditions: 


1. The only renormalizable operators that are present at the electroweak scale are those 
of the Standard Model, and the perturbative 4D effective field theory (dimensionally 
reduced from an intrinsically higher-dimensional theory) is valid. (The possibility of 
doing without an “X-boson” was briefly mentioned in |^.) In other words, we work 
with an effective theory with top as the heaviest field that has not been integrated 
out and neglect any Lorentz violating operators resulting from time dependence of the 
extra dimensional geometry. This condition is motivated from our ignorance about the 
physics beyond the Standard Model and the fact that no viable baryogenesis scenario 
has yet been proposed with only the Standard Model fields and gravity. 


2. The leading baryon number-violating operator due to physics at the fundamental scale 
is a dimension-6 operator of the form 


M? 


qqql 


( 1 ) 


where Vg the wave function overlap, which is a time dependent function. The wave 
function overlap Vg can be time-dependent whenever fermion localization is controlled 
by a field, akin to the radion, which is initially displaced from the minimum of its 
effective potential. The work of can also be taken to be of this type, where the 
“radion” is always at the minimum of a time dependent effective potential, which 
is derived from time dependent hnite temperature effects. Note that this condition 
implies the assumption that 4D instanton/sphaleron induced baryon number violating 
operators are weaker than the dimension-6 operator, i.e. only perturbative physics 
is analyzed for dimension 4 operators. This is a reasonable assumption as long as 
the temperature is well below the electroweak scale (which is anyway required in low 
fundamental scale models) and the SU(2 )l coupling remains small (Appendix P). 


3. The 4D Planck scale Mpi has an adiabatic time dependence due to the change in the 
volume of extra dimensions, and the fundamental scale of gravitational physics is at 
Mf ~ 10^ GeV. The upper bound on this scale comes from the desire to ameliorate the 
hierarchy problem without SUSY, while the lower bound is obtained by the bounds on 
neutron-antineutron oscillations, which cannot be suppressed by a geometrical mech- 
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anism. Note that there is still a certain amount of hne tuning, of order 10^, with such 
a (relatively) large fundamental scale. 

4. There is no electroweak phase transition because of a low reheating temperature {i.e. 
T < 30GeV) i, i- Baryogenesis is also assumed to occur above the QCD phase 
transition temperature {i.e. T > 0.2 GeV). 

5. The effective 3D space is homogeneous and isotropic at all times when the temperature 
of the universe is T > 0.2 GeV. 


In this context, we analyze the possibility of baryogenesis and hnd specihc conditions under 
which it is feasible. 

Our setup contains all the necessary ingredients of baryogenesis : 


1. Baryon number violation (dimension-6 operator) 

2. C and CP violation (Standard Model GKM -|- dimension-6 operator) - the resulting 
quantity is much larger than the Jarlskog invariant. 

3. Nonequilibrium (spacetime expansion or other unspecihed means) 


The reason for the difficulty of this scenario will be ingredient 1, the baryon number violating 
operator being too weak in the low temperature setting, independently of ingredients 2 and 
3. Since the usual Jarlskog invariant of the SM is not the relevant quantity, given the 
presence of dimension-6 operators, the reason why this baryogenesis scenario is difficult has 
nothing to do with traditional wisdom. In particular, we hnd that dimension-6 operators are 
in most cases too weak to compete with baryon number-conserving operators which create 
entropy, for most sets of reactants. 

Naively, one hnds this result a little surprising because one would expect 


_ nB 6cpvBn 
Vb = — ~- 7^ 

s g*sT: 


( 2 ) 


can be easily engineered to be sufficiently large, where Jcpv is the dimensionless quantity 
associated with CP violation, B is the maximal branching ratio into the baryon number vio¬ 
lating channel hxed by our model, n is the number density of particles that can be converted 
into baryon number, is the number of degrees of freedom contributing to the entropy, 
and T* is the temperature at some hducial time after which the baryon number is conserved. 


5 





However, even though 5cpv can be naturally large in our scenario (not governed by the 
Jarlskog invariant), Bn/{g^sT^) cannot in general be made large because as the baryon 
number violating reactions are creating baryon asymmetry, the baryon number conserving 
counterparts are creating proportionately large amounts of entropy density. In other words, 
having the quark-lepton separation go to zero does not generate order 1 branching of the 
baryon number violating channels. Instead, such channel is always suppressed by a scale 
much larger than the electroweak scale. 

We do hnd an exception to our “no go” rule for a special situation in which the universe 
starts out dominated by the b and c quarks and gluons (which might happen through the 
peculiarities of reheating after inflation), and the lighter quark c is out of equilibrium while 
the heavier quark b is in equilibrium, for example due to the anomalously large expansion rate 
and the anomalously small density of c quarks. In this exceptional scenario, there are two key 
ingredients compensating the suppression of the dimension 6 operator: 1. Due to the small 
iWfeP, the b quark interactions with out of equilibrium c quarks do not generate signihcant 
entropy; 2. The integration time for baryogenesis can be made long if the expansion rate 
is small, due to the nonstandard cosmology of extra dimensions. However, this remote 
possibility, which requires additional assumptions about special initial conditions and/or 
more new physics, is unlikely to be realizable, and we do not have a detailed model to 
demonstrate its full viability. 

Throughout this paper, we do not assume thermal equilibrium initial conditions, except 
where explicitly noted. Hence, although the short distance operators are hxed, the nonzero 
initial particle densities that exist in the universe can be a priori arbitrary combinations of 
species. Thus, when we consider bounds, we explicitly check all likely combinations of species 
without any bias to their initial number distributions, except when otherwise noted. As far 
as the out-of-equilibrium-condition due to the expansion of the universe is concerned, we 
assume that somehow the Planck scale can be adjusted to dial into the appropriate expansion 
rate. Of course, in a more completely specihed model, this may turn out to be problematic, 
due to the bounds on moduli overproduction and other light moduli problems. However, our 
results are robust with respect to how the moduli problems and initial conditions problems 
can be solved. In fact most of our results do not depend on how the particles are out of 
equilibrium. This broad insensitivity with respect to cosmological details is one of the key 
features to our work. 
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The order of our presentation will be as follows. We begin with a discussion of the class 
of braneworlds and the energy scales and couplings that we are interested in. We then 
present our constraint on this scenario for baryogenesis with dimension-6 operators in the 
low temperature regime. Finally, we summarize and conclude. 


II. BRANEWORLD SETUP AND SCALES 


The simplest way to realize the proposal of geometrical proton decay suppression is for 
the “brane” on which we live to be a nontrivial scalar held prohle in an extra dimension . 


On a smooth manifold, chiral fermions are localized in the presence of a scalar held prohle 
with Yukawa couplings to the fermion helds, resulting in a position-dependent ehective mass 
term. The simplest ansatz is to take a linearly varying scalar held, which results in chiral 
fermion zero modes with a Gaussian wavefunction, localized about the point where the 


ehective fermion mass vanishes 0. On an orbifold one may impose special boundary 
conditions, which themselves ensure that the lowest mode is chiral, then with a scalar held 
which is forced to take up a kink conhguration one obtains a fermion wavefunction which 
varies approximately as cosh“^(/ca; 5 ). 

By coupling the leptons and quarks diherently to diherent scalar helds, one may separate 
the zero mode wave functions such that the wave function overlap factor Vg in Eq. (|I]) 
obtains an exponentially suppressed value. As discussed in pl[|, it is not easy to induce large 
displacements of the points of localization away from special points such as the orbifold hxed 
points (whose boundary conditions force the zero of the scalar held) without introducing 
arbitrary scale parameters. In this paper, we will not assume any particular mechanism to 
make these localization points (which we will sometimes loosely call “branes”) dynamical, 
but rather suppose that large dynamical displacements of the localization points (branes) 
from their hnal value are possible. Hence, we treat Vg as merely a time-dependent function 
that is controlled by an exponential of the form 


Vg ~ e 


-CMfL 


( 3 ) 


where C is a constant coefficient of order 1 and L is a time-dependent function giving the 
distance between the quark and lepton localizations. There may be additional contributions 
to Vg varying as but these will be subleading unless jj? <C Mj. We will see below 
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FIG. 1: Schematic setup of our braneworld. 


that C, actually has to be greater than or equal to about 3/2 in our scenario. 

The baryogenesis scenario is then as follows: at some early time in the evolution of the 
Universe, M^L ~ C^(l) and the baryon number violating operator is suppressed only by 
the fundamental scale. By the time of nucleosynthesis, C^fL will have settled down to 
today’s value of about 50, consistent with proton decay constraints (see Appendix 0. For 
the purposes of our paper, we will not need the specific details of how the quark-lepton 
separation is stabilized at its final value. We will show that the short-distance physics 
prevents sufficient baryogenesis in the minimal model that we consider, except in the special 
case already mentioned. 

Since the zero mode quark and lepton wave functions are separated in the extra dimension 
and since the quarks and the leptons interact through the SU (2) l x U{1)y gauge group, these 
gauge fields must also propagate in the extra dimension, and the 4D gauge theory must come 
from dimensional reduction. Similarly, since everything interacts with gravity, 4D gravity 
also arises from dimensional reduction. We now derive bounds on the fundamental scale 
as well as the size of the extra dimensions by demanding perturbativity of gauge helds 
propagating in the extra dimensions and that Einstein gravity be recovered 0. 

From gravitational dimensional reduction, given that there are a total of n extra di¬ 
mensions and one of them is asymmetric with length the others having lengths i?*, 

i = 1,..., n — 1, we have 

/ur2 

/Vf2 ^_ l'4'l 

where ts denotes the present time and Mpi is the reduced Planck mass Mpi/\f^. There is 
also a perturbativity constraint on the Yang-Mills coupling constant. In 4 -|- p dimensions 






(for example for a (3 +p)-brane), the Yang-Mills coupling is 

2 ^ ^4+p 

y^+p 

where A 4 +P is a dimensionless number. By dimensional reduction to 4D, we have the 4- 
dimensional dimensionless Yang-Mills coupling g 4 ^ related to A 4 +P by 

= (5) 

i 

If gA{t) is to remain perturbative at time ti < t^, we must have 



Given that 5 f|(t 3 )/( 47 r) 1/31 for the EW theory [a = ~ 1/137 = Q;vysin^ 6 *w — 

0.225aw], we require 

(n(f||))< 31. 

If we have a scenario in which 

Rbit) = const. 


this is not an issue, but in general, this is constraining. 

Now, let us consider absolute perturbativity of the hve dimensional theory. Since the 
cutoff is set at Mj, perturbativity requires 

< 1 
dvr 


Using Eq. (^, we obtain 

% = ^{Rdh)M,)\[{Ri(,h)M,) < 1 

i 


( 6 ) 


for 5D perturbativity. However, this requires the volume enhancement factor to be small. 
In the end, if we want to minimize this effect, we must have only a D4-brane (he. p = 1). 
Note that our perturbativity constraint is for calculational feasibility rather than necessarily 
for fundamental physical reasons. 

Coming back to Eq. (^), if we impose that in which the gauge helds propagate is 

the largest extra dimension 


Ri < R^{h) ^ 


1 

awMf 
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where the approximate equality is due to Eq.(|^), one sees that 


Mf > ^ 1.5 X lO^^GeV. 


To restate this, if the gauge helds come from dimensional reduction of the largest compact 
dimension, perturbativity requires 


Mf > 10^=^ GeV. 


Hence, for this scenario to have Mf < 10^ GeV, there must be compact dimensions larger 
than the dimension in which the gauge helds propagate. 

Then we take, for simplicity, 5 large extra dimensions of length R' in which the gauge helds 
do not propagate, and one comparatively smaller extra dimension of length R^it^) = 
in which (electroweak) gauge helds do propagate. Since in that case, 

2 _ awM^i 
f {R'Mff 


we can choose 

2 X 10^ /lO^GeVV''® ^ 1 /lO^GeV]^^^ 

Furthermore, the experimental limit from neutron-antineutron oscillations forces Mf > 
10^ GeV PI, therefore R' must be smaller than about 2GeV“^. In the spirit of the orig¬ 
inal motivation for large extra dimensions, namely in order to ameliorate the hierarchy 
problem, we take the smallest consistent value of Mf and set it to Mf ^ 10® GeV. 

What is important to keep in mind from this section is Eq. (!)• Since oiw{G) = ~ 

1/31, the maximum value that V{G)Mf can take is 31. Hence, proton decay constraint 
requires C, > 3/2 as claimed below Eq. (^. Furthermore, to “unsuppress” the dimension-6 
operator, LMf 1 is required during baryogenesis. This means that if the distance between 
the branes L is equal to the time dependent size of the extra dimension i? 5 , then aw ~ 1 
during baryogenesis. However, as shown in Fig. Q, if the distance between the branes can 
be adjusted independently of the size of the extra dimension, this need not be true. Finally, 
if we require that no K-K modes of the SU{2)l x U{1)y gauge group contribute to the 
baryogenesis dynamics, we must have > 100. This leads to the requirement that 
awit) > 10“®. We will henceforth take aw to be a free (time-dependent) parameter, with 
a value 

10-® < aw{t) < 1. (7) 
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With a setup in which proton decay and other 5-violating processes are suppressed by 
the quark-lepton separation, the next question to ask is which processes dominated at the 
epoch when L was smaller {L~^ ~ Mf). Proton decay due to massive fermion exchange, 
which appears in the effective theory by operators suppressed as with fc ~ 1, is 

the most “difficult” to suppress at the present epoch, since one requires MfL > 50 at the 
present time up to factors of order 1, and R is bounded from above by if SU(2)w is 
to remain perturbative in the bulk. Decay mediated by local interactions in the d = 5 
effective field theory, requiring the localized SM fermions to propagate through the bulk 
with a suppression of might also be problematic if /i^ is small relative to M|, as it 

must be for a credible field theory model of localization. Here we would require > 50 
today. Nonperturbative effects involving the extra dimensions should also be considered, 
for example virtual black holes or wormholes connecting the two branes. However, given 
very reasonable assumptions for the topology of the corresponding instantons, their action 
is large enough for the effects to be negligible [Q. All other i?-dependent effects turn out 
generically to be smaller, including the weak instantons which do not cause proton decay but 
might be important in baryogenesis (Appendix B- Then, assuming that both suppression 
factors saturate the current bound, as R shrinks from its present value the effective operators 
suppressed by grow much more quickly than those varying as g-const.M/H^ time 

of baryogenesis, for R comparable to Mf, the dominant processes (all other things being 
equal) will be those mediated by the SM fermions propagating through the bulk, for which 
the suppression is ~ e ^ ~ 1. 

The goal now is to see if without an electroweak phase transition, we can create enough 
baryon asymmetry with the dimension-6 operator which is unsuppressed. 


III. BARYOGENESIS 


The dehning characteristic of our class of baryogenesis scenarios is the leading baryon 
number-violating dimension-6 operator, one example of which is written in the gauge eigen¬ 
state basis as 


AT = 


UoA 


anbr 


M] 


(.QDic-'iQlWinc-'iQlUce 


JjJm 


+ + cae 


im jl] 


( 8 ) 
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where we have suppressed color indices, are family-dependent coefficients, {i,j,/,m} 

are SU{2)l indices, summation of indices is implicit, and ci, C 2 and C 3 are constant real 
coefficients. The wave function overlap integral that gives rise to Vg will be taken to be 1 
during baryogenesis. Rotating the quarks and leptons into mass eigenstates, we have 


{Q% = Uf{Q% 


(9) 


(L“), = Wf{L% ( 10 ) 

which corresponds to the usual flavor mixing matrices Vckm = U}jJd and Vmns = W^Wd- 
This then gives 


AL 


\anbr 

Mf 


(Ql)ic-'(QtULl)jc-'(Q 


LJm 




( 11 ) 


where 


yanbr _ ya'n'b'r' jja' ajjn'njjr'r'y^b'b 


( 12 ) 


Note that due to the introduction of the dimension -6 operators, the global symmetry of 
the theory (in the limit that the Yukawa couplings vanish) is different from the usual low 
energy analysis of the Standard Model. Hence, we have more than one rephasing invariant 
CP-violating complex parameter. 

Now consider for example the interference between the tree level and the one-loop dia¬ 
grams shown in Fig. ^ assuming for simplicity that C 3 = Ci = 0 and C 2 = 1. If the tree level 
amplitude is denoted as Ali while the one loop amplitude is labeled as ^ 2 , we know that 
M.I can be written as 

Ml = RiAfi '22 (13) 

where Ri is purely real, while M 2 will in general be of the form 


M 2 = 



(14) 


where R 2 and I 2 are real functions. Note that I 2 will in general arise if we can cut the 
propagators in the loop (to make two tree level graphs) and put them on shell simultaneously. 
For the CP conjugate reaction Mi = Ri\fi 22 = but 

M 2 = {R 2 + ^ M *2 (15) 


where in particular, the imaginary part arising from the cut propagators is not the complex 
conjugate of that for the CP-conjugate reaction, since it derives from kinematics alone. 
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FIG. 2; Diagrams interfering for baryogenesis. 

In the Boltzmann equations for the baryon number density, the probability of interest 
will be 


D = = 2Re|{A<,A<* - MiMD] 

— 4 { Ri (/2 + 12))^ 


(16) 

( 17 ) 


where we have defined 


d= (is) 

and the {) represents spin averaging. Written out explicitly, this interference term is 


D = 


/ o 


{^ 1122 * ^ii 22 VisVtj}IsUj Re {Tr [^3 7 ^ 7 ^^! ^4 lulph Pr ]} + Pai’ity odd (19) 


where the external momenta Pi are labeled clockwise and constrained to satisfy Pi + P 2 ~ 
P 3 -P 4 = 0, Pr= {1 + 75 )/ 2 , and 


rO/3 _ 
^stij — 


l2x / 

(P3 - A;)“ 

\ {pa + kY 

1 1 

1 J (2>r)‘‘ 

_{P3 - kf -ml. _ 

1 

(M 

+ 

-It 

3 

k^-Mly\ 


( 20 ) 


Note that vanishes only if the intermediate quark states cannot go on shell. If the 
intermediate quark states can go on shell, the expression becomes 

d^qi dPq2 


1 % = -Stt' 


-{2TTY5^^\qi + q 2 - {pi +P 2 ))- 


Q.iq.2 


( 21 ) 


(27r)32E,, (27r)32E,, ^ ^ {p, - q.f - ' 

Now, let us make some estimates. Consider the standard parameterization of Vckm- 
Suppose we choose the angles as 6*1 = 3.366, 62 = 0.0307, 6*3 = 0.0067, and a phase (5 ckm = 
1.4633 roughly consistent with experimentally known bounds. We specify in part a Yukawa 
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FIG. 3: A typical distribution of CP violating coefficients from a random texture. 


texture by choosing the rotation matrix 




^ 1 0.5 0.9^ 


Uu = exp 

i 

0.5 1 0.6 




^0.9 0.6 -2y 



and suppose Ui = Wi. With a array of random (real) numbers of order 1 for the we 

hnd the effective CP violating imaginary coefficient to be distributed as shown in Fig. 
where we have summed over the indices i, j as if rest of the transition probability did not 
depend on i and j. As expected, the quantity 6 (analogous to the Jarlskog invariant in the 
presence of the dimension-6 operators) is around 0.1, much larger than the Jarlskog invariant 
of ~ 10“®. We expect similar results if the are taken complex with random phases. 

As for the we have when the Mandelstam invariant s 

Isfij ~ ( 22 ) 

while for the top quark, 

~ 30. (23) 


Note that Mw can be time dependent since 


Mw ~ \/2Tiawv 


(24) 


and aw can be time dependent (anywhere between 0 < aw < !)• Although v may be 
time-dependent, we will take it to be approximately constant [pl|] . 


14 









Now, let us write down the Boltzmann equation for our system of quarks and leptons. 
Denoting the phase space density of particles in an FRW cosmology with scale factor a as 
fx, we have 


dfx a 1^2 Qf^ 


dt a Ex dEx 


= E J\-Ui>\Mi\^-2fxhRe{MiM;} + 


fcfd 






fcfd 
req peq 




feq ft 

jc Jd 

~fxfb\-M. 3 \‘^ ~ 2/a;/6 Re + 

f^fj 


fifj 

fcg yeq 


fx^fb^lMsl 


-2 


r/r 


rjn^ReiMsMl} 


where 


2Ex. 


^ereactants 


E Pr 

rereactants 


(25) 


(26) 


is an operator generating an appropriate momentum integral (with the measure dlf^ = 
dfpn/[2En{27iY]), M-i is the amplitude for the tree level process xb —> cd, and M .2 is 
associated with the one loop correction with intermediate states ij. The amplitude is 
associated with the tree level amplitude of xb —>• ij, where ij are the intermediate states in 
the calculation of Af 2 , and the one loop correction to Afs involving the intermediate states 
cd (appearing in Afi) is denoted as M. 4 . Here, we have neglected the Pauli blocking factors 
for the Boltzmann collision integral. Although this is easily generalizable to include the 
Pauli blocking factors and although it is not justihed when calculating an accurate number 
density evolution, the neglect is justified when calculating an upper bound on the baryon 
numbers produced. Hence, for clarity of physics we shall drop them in this paper. Because 
only the amplitudes with the fermion propagators cut in Ai 2 and A 44 contribute, one can 
show that within the collision integral sum that RejAfiAl^} = — RejAdsAl^}. Note that 
the collision integral vanishes in equilibrium. (A related work on this kind of Boltzmann 
equation may be found in Ref. |2^.) 

Thus one arrives at the baryon number density evolution equation 


1 dijiBof) 
dt 



fifj + fjf-j 

f^fr 


fcfd + fcfd 


fcq fcq 

Jc If 


(27) 

where ba is the baryon number associated with reactant a and the prime on the integral 
denotes that all the momenta are to be integrated over. Note that for each baryon number 
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associated with reactant a to contribute to the baryon asymmetry, the intermediate and the 
hnal state must be out of equilibrium and not the initial reactant. This is expected from 
unitarity and CPT which implies that 

(28) 

i j 

where j sums over all hnal states: when the probability is summed over all hnal states, there 
cannot be any CP-violating enhancement or subtraction of the probability that survives. 
Also, note that although naively one might think that the RHS of Eq. (|^) would vanish, 
since the quantity inside the square bracket is odd under the exchange ij •«-*• cd and the 
summation is over all h j, c, d, there is no cancellation because Re C cd)} 

is also odd under the exchange ij ^ cd. Finally, note that although we have only written 
down the scattering processes explicitly, the decay and inverse decay reactions have an 
analogous form. 

Now, instead of accounting for the possible contributions to the baryon asymmetry from 
all the reactants together, let us put bounds on the the maximum contribution that arises 
from each of the reactants. As explained in Appendix we can derive from Eq. (^) 

2rj 


where 


r 

riB < — dtTBvit)a^nf, 

CXj \J 


11 = awY,1 1 “ ~ 0.02 


f2 


(29) 


(30) 




~ — Tinu and rriinit are temperature and mass scales, respectively, 

of the initial state (or hnal state) reactants, and rudj are the masses of the particles in 
the loop, Tbv is the tree-level baryon number violating reaction rate which can be read off 
from Table |. In the Table the function W dehned as 




(31) 


+ A{T/mQaY 

is a time dilation function arising from the thermal averaging and F{mQ,mB,TQ,TL) is an 
function of order 1 also arising from thermal averaging {Ki is a modihed Bessel function of 
the second kind). Using these and other results presented in Appendix y, we can place an 
upper bound on the baryon asymmetry generated by any single species out of equilibrium 
as 


Hb 


< 27) 


rBv(^) 

rtot(^) 


(32) 
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Fbv 

reaction 

magnitude 

AiB + L) 

FbV dec(l + e) 

Ql ^ L\QlQi 


-2 

Fbv dec(l — e) 



2 

Fbv scatt(l + e) 

L^Ql ^ Q^^Qi 

nL,Qp {rriK +T^ +Tl ) 

--(l + e)F(mQa,m^6,TQ,rL) 

/ “ 

-2 

Fbv scatt(l — e) 

mi ^ QlQi 

+t 2+r2) 

--(1 e)F{mga,mLb,Tg,TL) 

f “ 

2 


TABLE I: The dominant reaction rates. The upper indices on leptons L and quarks Q label family 
number while the subscript indices label denote up-type or down-type {e.g. is a neutrino). All 
other quantum numbers are implicit. The small number e indicates the splitting due to the CP- 
violating phases. Only the quark decay rates are displayed, as they will be the most relevant 
quantities. The coupling A of the dimension-6 operator is assumed to be of order 1. 


where the maximum is taken at any time during baryogenesis. Note that this bound is 
independent of cosmology (without specifying the time dependence of the scale factor a) as 
long as the usual Boltzmann equations are valid. 

Let us hrst consider the quarks. Instead of writing Ftot as the sum of rates for all the 
channels, it is more illuminating to consider separately the cases when Ftot is dominated by 
either the baryon number-conserving annihilation rate or the decay rate. The 5-conserving 
annihilation rate is dominated by a t-channel quark annihilation into gluons, with the pos¬ 
sible exception of the top quark for which the rate can be dominated by a weak scattering. 
Explicitly, we use 

2 

hann ~ nq F2(mQ, Tq) ( 33 ) 

Tt-ann ~ Tann ( f —^ |f^CKMp = Tann/^i ( 34 ) 

V Mw J \Ois J 


where |VckmP is the appropriate quark mixing matrix element squared, and F 2 is a function 
encoding the effect of thermal averaging, similar to F. The 5-conserving decay rate is 
dominated by a weak decay rate, which is enhanced for the top quark due to the longitudinal 
W component going on shell. We use 


weak dec 


dec 


atv f rnq 


1671 \Mw 


mQ\VcKM\^W {mga^Tq) 


TT 


weak dec" 


Qiw \ 


Ml 
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= F 


weak dec 


/F, 


(35) 

(36) 
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ratio # 

FBv/Ftot 

1 s)max 

1 

FbV dec/Fann 

2rj \mQ(m'^Q+T^)W(mQ,TQ)'\ t.. \ 

^ nQF2{mQ,TQ) v 

^ J 1 J max 

2 

rBV scatt /rann 

V (™Q+T3)(™/,+"1q+T3) 

nf^(t)F{m.Q,m.f^,TQ)-\ 7 ,, 

00 

p 

nQ(t)F2{mQ,TQ} 

3 

TbY dec/^weak dec 

2ri t Mw 74 1 t TP\ 

4 

^BV scatt/r^weak dec 

327777 7 Mw 74 +™Q+'^Q+^/2 ^ 

g'^ ^ Mf ' m^IVcKMk 

n/,jF{777Q,m^2TQT/2)] 7 .. 17 7 


TABLE II: nsjs for quarks. For the top quark, each of these ratios must be multiplied by Et or 
Ft- However, these extra factors do not change our conclusions. 


Explicitly, Et ~ ( igocev ) ''^^ile Ft ~ 0.05 ( igQQev ) ’ which implies that the weak 

decay rate is not much enhanced for the top quark through the Goldstone emission. Given 
these rates, we list the leading bound hb/s [i.e. the upper bound on the baryon fraction 
relevant to the case that the given reactions are the dominant contributions) in Table 
The lepton sector is similar except with Ug —^ Oem and no top quark enhancements. In addi¬ 
tion, the leading number density changing channel for the neutrinos may be coannihilations 
instead of self-annihilations. From the last section, we have Mf > 100 TeV which implies 

(37) 


Finally, we need to constrain the number density of the particles. Since the elastic scattering 
rates are at least as large as the total inelastic scattering rates, to absolutely forbid any phase 
transitions, we restrict the number density of particles to be smaller than the thermal density 
at T* = 30 GeV, i.e. 

ngj < 0.2Ti> (38) 


if relativistic or 


< 0.13(mo,ir.)’''"e-'"«’‘'’'- 


(39) 


if nonrelativistic 32 


Now, consider the case when the branching ratio is approximated by ratio 3 in Table ||. 
For the top quarks out of equilibrium, we have 


ratio 3 = 5 X 10 
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180 GeV 
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(40) 


18 






























while for all other reactants, we have 


ratio 3 < 6 X 10 


-17 


\ 2 


(41) 


VI^ckmIV VlSOGeV. 

where we have used the bottom quark mass of about 4.3 GeV. For the top quark, because 
the upper bound here is smaller than what is needed for sufficient baryogenesis and since 
the same numerator appears in ratio 1, the 5-violating decay channel is ruled out for the 
top quark. For the bottom quark, since the leading decay channel has \Vch\^ ~ 1.6 x 10“^, 
Eq. ( 1^) also rules out sufficient baryogenesis. All other channels also fail the tests under 
ratios 1 and 3 which means that the S-violating decay channel cannot help in generating 
sufficient baryon asymmetry. 

Gonsider now ratios 2 and 4. Ratio 4 requires for the top quark that 

riB 


< 4 X 10 


-15 


(42) 


if we assume that F ~ hF ~ 1. Hence, ratios 3 and 4 completely rule out the top quark as 
a signihcant contributor to the baryon asymmetry. 

Let us now consider the scattering of the out-of-equilibrium particle / with every other 
particle labeled as /2 in view of ratio 2. It can easily be seen that ratio 2 is not very 
constraining, as it can in principle be made arbitrarily large. Hence, we can consider a more 
stringent version of the entropy bound corresponding to ratio 2, Eq. which we rewrite 
as 


hb ^ + m\ +T] + T]^) 1 ftf 


Mjs 


JtB 


dt Tf^)nfnf^ 


(43) 


where s is the entropy density, Ti is the temperature of particle species i, Tj is the charac¬ 
teristic temperature of species i inside the integral, gi is a kinematic function of order 1, 
Ib marks the beginning time of baryogenesis, and tf marks the end time of baryogenesis. 
Instead of using the maximum contribution to the entropy, it is sufficient to consider any 
one contribution since we are placing an upper bound. Since we are guaranteed that the 
reaction changing the number density of / is out of equilibrium, and there exists a contri¬ 
bution to the entropy from the weak interaction (//2 ^ hnal products) with the matching 
integral of the form dta^UfUf^ (since / is necessarily out of equilibrium), we can write 
using Eq. 

s > 


«wl^CKMp(^/ + + Tf -F 1 rtf 




JtB 


[ ^ dta^g2{mf,mf^,Tf,Tf^)nfnf^ (44) 
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where g 2 is another order 1 dimensionless fnnction. It is important to note that we have 
assnmed the kinematic viability of the weak interaction in Eq. Combining Eq. 

and Eq. (^), we have 
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27] 


- al^lVc 


M, 
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KM 
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5 X 10 


-13 
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180 


V y 

G^J 


(45) 


I VckmI 

independently of aw- Hence, with |VckmP < 10“^, snfficient baryon asymmetry generation 
seems feasible. A case in which Eq. (^5]) applies is when ont of eqnilibrinm b qnarks are 
interacting with the c qnarks at an effective temperatnre of 30 GeV. In particular, if the 
baryon number conserving weak interaction channel is dominated by be —>■ cs which gives 
the iHcfep ~ 1.6 X 10“^ suppression, sufficient baryogenesis in this scenario seems possible. 
We will return to this possibility later and constrain it further. 

Another possibility to evade Eq. (|45|) is to have the weak interaction be kinematically 
suppressed. In the case of out of equilibrium c quark, the reaction cUr —*■ st has the 
kinematic suppression necessary to evade Eq. (^5[) in the limit of zero temperature since 
rrij. ^ 1.7 GeV while rric ~ 1.5 GeV. In Table B we list all the weak interactions that are 
kinematically suppressed such that sufficient baryogenesis is possible with respect to the 
constraints thus far. 

According to Table B some weak interactions are kinematically suppressed only below 
the temperature of the QGD phase transition, which is assumed to occur at around T 0.2 
GeV These interactions hence really are not suppressed in the temperature range 

of interest and the baryogenesis bound is not evaded for these reactants. Hence, these 
are labeled irrelevant in Table H- Only interactions involving Ur then remain as viable 
candidates for the kinematically suppressed baryon number conserving channel, if we assume 
that its mixing with or is less than 10“^. If the neutrino sector is mixed, all the bounds 
return and there are no candidates that can generate sufficient baryon asymmetry . 

In both cases where Eq. (|4^) allows sufficient baryogenesis, the baryon number-conserving 


weak interactions between the reactants contribute too little to the entropy to rule out 
baryogenesis. In such situations, let us consider whether the self-annihilation interaction 
contribution to the entropy is enough to dilute away the baryon asymmetry contribution. 
Start with Eq. 


and again approximate ~ 1. However, instead of Eq. (M), consider 

^ ^ f‘< (4g) 

a'^ Ju ■' O'* Ar ■' 
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out of equilibrium particle 

weak interactions 

Threshold Temp. (GeV) 

Relevant 

e.g. BV 

c 

CVr sf 

0.4 

Y 

cvr dd 

s 

SUr — > UT 

1.5 

Y 

SVr du 

d 

dUr UT 

1.7 

Y 

dVr du 


du^ —> up 

0.1 

N 


u 

uVt — > df 

1.7 

Y 

uVt — > dd 

55 

— > dfi 

0.1 

N 


55 

uue — > de 

0.006 

N 


e 

ue dve 

0.005 

N 


Vr 

u-rC — > dr 

0.21 

N 


55 

V-rS — > UT 

1.505 

Y 

UrS — > du 

55 

Urd — > UT 

1.7 

Y 

U-rd —> du 

55 

U-rU —> dT 

1.7 

Y 

UrU dd 


u^d —> up 

0.1 

N 


55 

—> dp 

0.1 

N 


Ve 

VeU —> de 

0.005 

N 



TABLE III: The set of reactions evading the simplest bound. These reactions are suppressed when 
the temperature falls below the threshold temperature. 

which corresponds to the entropy bound arising from self-annihilation reaction contributions, 
where the contribution to the entropy from the species /2 only begins once it goes out 
of equilibrium, at time tc, while the / self-annihilation contribution persists throughout 
the entire duration [ts,t/] of baryogenesis. Here, the thermal averaged quantities (such as 
(aAfv)) have their temperature hxed at some value, characteristic of the integration interval. 
The advantage of considering the bound arising from self-annihilation reactions is that they 
change the number density and are always possible at tree level if there is a neutrino lighter 
than (or degenerate with) Ur- We shall assume that the self-annihilation channel is open 
Suppose first that the term dominates over the second term. Then, using the 
elementary inequality 

n% +n}> 2nfnf^ ( 47 ) 
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we can write 


where we have dehned 


Hb 
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s 2M‘j{aAf2v) 


(aAf v) 

1 + y^ + e 

WAfV) 


and 


M] 

_ /,*• dt a^n% 
/,'/ dt 


( 48 ) 


(49) 


(50) 


which measures how much the /2 particles can contribute to baryogenesis through interaction 
with out of equilibrium particle / without contributing to the entropy. This can be seen from 
noting that Uf^ is the equilibrium distribution during t G [tB,tc] by dehnition. Producing 
baryon asymmetry for a long time without producing any entropy (i.e. letting tc approach 
tf) is the best that can be done to circumvent the suppression of the dimension 6 operator. 
For sufficient baryogenesis {jib/s ~ 10“^°), we thus require 


e> 


i^AfV) 


where 


5 = 
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MRcTAfV) 


X 10 


10 


If the rij dominates in Eq. (|4^) , one can similarly derive Eq. 
is true in general. Now, noting that one can write 




and using Eq. (|47|) , one can write 


(2 — ( 5 ) f dta^nfnf^< f dta^nj^. 

JIr JtR 


It follows that 


,3 (1 + 7 ) 

/ dta nfUf^ < 

Itg 


- -— / dt a?n^f 

(2 - 4) Jt, 


(51) 

(52) 

thereby showing that it 

(53) 

(54) 

(55) 


which is useful if ^ > 1 and 5 < 1, since as we mentioned before, nf^ is the equilibrium 
distribution during the time interval [tB, tc]- Hence, we shall use Eq. (^) with parameterized 
cosmology to put a bound on cosmology from the condition of sufficient baryogenesis. 

We parameterize the unknown cosmology as 


a (X t"'. 


(56) 
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(57) 


and determine tc by the ont of eqnilibrinm condition 


a(tc) 

a{tc) 


= max [j {aAf 2 v)nf^ (4), T/J 


(58) 


where j is an 0{1) constant and Fis the decay rate of f 2 - Becanse the ratio of the decay 
rate Fto the expansion rate increases as a fnnction of time, when the decay becomes the 
dominant chemical potential changing reaction, it can lead to an ont-of-eqnilibrinm nnmber 
density of the /2 only if the hnal states of the decay are not in eqnilibrinm. Hence, the 
maximnm above applies only if the decay prodncts are ont of eqnilibrinm, and otherwise, 
only the self-annihilation shonld be taken into acconnt. Fnrthermore, note that any choice 
of tc that we make in Eq. (^8D corresponds to a different choice of cosmology becanse 4 
sets the scale of expansion rate of the nniverse. The larger tc is, the slower is the expansion 
rate, and hence the slower the dilntion rate of any accnmnlation of baryon asymmetry. As 
far as interpreting the parameter u is concerned, it is shown in Appendix 0 that z/ < 1 for 
most physical systems (even with external sonrces snch as the inflaton) conhned to 3 spatial 
dimensions, nnless there exists a cooling mechanism different from the expansion of the 3 
spatial dimensions. However, with the presence of extra dimensions, the energy exchange 
with the bnlk exotics may indnce a more rapid cooling allowing for u > 1. If {o'Af 2 'v)nf^(tc) 
determines the ont of eqnilibrinm time of / 2 , we hnd 
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s jcMj{aAf 2 v) Tf (2 - 6) z/tt^ Jt,/Tb 


1 /y2 ( ™'/2 . y ') 


A'2(=f) 


(59) 


where c = 27r^5f*5'/45, is the nnmber of degrees of freedom contribnting to the entropy 
at time tc, and the factor of C 2 acconnts for (the sqnare of) the degrees of freedom in the /2 
nnmber density {K 2 is a modihed Bessel fnnction of the second kind). We will assnme that 
acconnts for at least the 16 degrees of freedom of the glnons in the thermal bath. As 
for C 2 , if /2 for example is a qnark, then it will be 9 for the 3 color degrees of freedom (the 
2 helicity states were already acconnted for, which means that if /2 were Weyl nentrinos, 
this shonld be divided by 4). If Fdetermines the ont of eqnilibrinm time of / 2 , then this 
qnantity should be multiplied by 


{aAf2v)m\TK2{^) 

7r^r/2 


(60) 


23 






For the relevant reactions in Table |T|, a Vr and a qnark are always paired. For the qnarks 
ont of eqnilibrinm, rrif^ = rriy^ ^ 0 and m/ = mq\ 






JK + T, 




- 1 
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2r]{ml + ml 


+ f^ + fl){ml + T, 
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alMj 


-1, 

(61) 

X 10^°. 

(62) 


These eqnations tell ns that for all effective temperatnres (between 1.7 and 0.2 GeV) and 
qnark masses (less than 2 GeV) relevant to onr scenario, we have the conditions 5 -C 1 and 
3> 1, which allow nontrivial bonnds to come from Eq. (|^. From Eq. (|5^), it is easy to 
show that there are no valne of i/ > 0 and n > 0 snch that snfficient baryon asymmetry is 
generated. Hence, there are no scenarios in which there is snfficient baryogenesis throngh 
ont of eqnilibrinm qnarks interacting with Ur- 

For the case of ont of eqnilibrinm and q in eqnilibrinm, mf^ = mg and Tnj = k. 0 
and we have a similar resnlt; 
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(63) 

(64) 


Again, it is easy to check that for all relevant effective temperatures (between 1.7 and 
0.2 GeV) and quark masses (less than 2 GeV), we have the conditions 5 -C 1 and ^ 3> 1 which 
allow nontrivial bounds to come from Eq. (|^). As before, Eq. (|59| ) rules out the possibility 
of sufficient baryogenesis for n > 0 and n > 0. Hence, even neglecting the mixings, it seems 
nearly impossible to have sufficient baryogenesis with Ur through the dimension-6 operator. 

Let us now return to the case where the baryon number conserving weak interaction 
channel is dominated by be —> cs which was seen to have sufficient suppression for us to 
further consider the possibility of successful baryogenesis. In applying the reasoning that 
led to Eq. (^^, we can improve the bounds compared to just Eq. (|5^) by considering that 
the entropy also has a lower bound at temperature of 


s(Tc)a^iTc) > ng{tB)a^{tB) + 


(65) 


where Ug corresponds to the density of gluons and by putting R inside the time integral and 
tying all the temperature scales together. The bound then becomes (similarly to before) the 
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FIG. 4: A conservative upper bound for the baryon asymmetry generated by out of equilibrium 
b quarks interacting with c quarks. The parameters n and v parameterize various possible non¬ 
standard cosmologies. A value n > 1 corresponds to inflating cosmologies and in the absence of 
abnormal cooling v <1. This scenario does not generate enough baryon asymmetry. 


smaller of either 
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( 66 ) 


s uM‘j{aAf 2 v) Mj 


or Eq. (^). The large enhancement for the b-c interaction scenario in both Eq. (15^) and 
Eq. (|6^) comes from the long time tc during which baryogenesis occurs while /2 is in equilib¬ 
rium. This leads to an accumulation of any net gain in baryon asymmetry over the entropy. 
More explicitly, tc varies as 1/K 2 {mf^/Tc) (or the equivalent factor from Eq. (|60|) when ap¬ 
propriate) which comes from the exponentially large ratio of the short distance physics scale 
to the expansion rate scale when /2 goes out of equilibrium. The expansion rate must be 
exponentially small when /2 goes out of equilibrium (if at all) because the self-annihilation 
reaction must maintain equilibrium even in the nonrelativistic regime in which the number 
densities are exponentially suppressed. The resulting bound with mf = mb = 4.3 GeV and 
mf^ = me = 1.5 GeV is shown in Fig. As one can see, even with such an enhancement, 
there is insufficient baryogenesis (or marginally enough) in most cases. Even for the marginal 
cases, we do not expect the scenario to be realizable given the conservative nature of the 
bound. Other scenarios with small \Vckm\‘^ weak interactions while the heavier particle is 
out of equilibrium {i.e. out of equilibrium b interacting with u, as well as out of equilibrium 
c interacting with d) fail to generate sufficient baryogenesis for similar reasons. 

If we do not assume thermal equilibrium initial conditions (before relevant species go 
out of equilibrium), then as before we may consider the b being in equilibrium while the 
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c quark is out of equilibrium. According to our explanation of the enhancement for this 
scenario before, the larger the mass of the particle that stays in equilibrium until temperature 
drops to around 0.2 GeV, the smaller the exponentially suppressed density of nonrelativistic 
particles, and then as long as the self-annihilation determines the out of equilibrium time 
tc, the longer is the time tc over which baryon number is generated, leading to a larger 
accumulation of baryon asymmetry. Since mh > rUc in this scenario, we thus expect larger 
accumulation of baryon asymmetry over the entropy in this scenario. In fact, the self- 
annihilation reaction becomes small by the temperature T = 0.2 GeV to the extent that 
it no longer controls the out of equilibrium time (instead it is determined by the 6-decay). 
The competing entropy generation is suppressed because of \Vcb\'^ ~ 10“^ for the weak 
interactions, and other reactions generating entropy are insufficient as well in diluting the 
baryon asymmetry. 

For completeness, one should note that since irtf < in this scenario there is a contri- 

). This comes from the fact 


bution to the entropy not accounted for by Eq. (|59D and Eq. 
that 6 undergoes out-of-equilibrium decay to c even though 6 has a chemical equilibrium 
distribution. That is because c by dehnition is out of equilibrium and thus cannot give the 
inverse decay reaction necessary for the 6 decay to conserve entropy. However, we can check 
explicitly this contribution to the entropy contribution is negligible, even though it is pro¬ 
portional to the integration time because of the smallness of the decay rate. Also, since 
the c is out of equilibrium, the reaction of 6-6 annihilation into c-c will also contribute to the 
creation of new entropy. However, one can check that his contribution is also insufficient to 
dilute away the baryon asymmetry generated. The numerical bounds shown for this scenario 
in Fig. § indicate that this scenario is viable at the level of our analysis. 

The most obvious difficulty with this scenario is that another mechanism is required to 
set the initial density of 6 quarks to be much larger than that of c, as well as to maintain 
that hierarchy. For example, for temperatures far above rUc, the gluons which keep 6 quarks 
in equilibrium will tend to produce a large number of c quarks that will be close to chemical 
equilibrium. Since a high initial temperature is crucial, some mechanism besides just the 
expansion of the universe is required to keep c out of equilibrium while is 6 in equilibrium. 
Even if the c densities could somehow be kept small, since the self-annihilation cross sections 
for 6 and c are approximately the same, if the cosmological expansion is responsible for the 
out-of-equilibrium condition, then it must delicately tune itself to be smaller than the 6 
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FIG. 5: A conservative upper bound for the baryon asymmetry generated by out of equilibrium 
c quarks interacting with b quarks. The parameters n and i' are as before. For nonstandard 
cosmologies, sufficient baryogenesis may be possible. However, as explained in the text, the out- 
of-equilibrium condition seems difficult to maintain. 


scattering rate while larger than that of the c’s. However, our “model-independent” analysis 
does not rule out this scenario, and it may be viable in a complicated and fortuitous setting. 

As for the scenario in which b quarks are in equilibrium while u quarks are out of equilib¬ 
rium, it will be even more difficult (if not impossible) to keep u quarks out of equilibrium for 
the required “exponentially long” duration tc, because the gluons which are supposed to be 
in equilibrium have one of their strongest equilibration reactions through u-u creation . 


For the scenario with other quarks (besides the top and the bottom) in equilibrium during 
baryogenesis, the exponential enhancement of the duration of baryogenesis is too weak to 
allow sufficient baryogenesis. 

Hence, remarkably, only one scenario, namely out of equilibrium c quarks interacting 
with equilibrium b quarks, may generate sufficient baryon asymmetry with the dimension 6 
operator if the out of equilibrium mechanism can be appropriately engineered. However, it 
would be surprising if such a realistic scenario can be constructed. 


IV. CONCLUSION 

We have studied a new scenario of baryogenesis involving dimension-6 baryon number- 
violating operators that would generically arise from integrating out ultraviolet degrees of 
freedom in the context of the Standard Model. The question to be answered in this scenario 
was whether a time-dependent dimensionless suppression coefficient for the dimension-6 op- 
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erator can allow sufficient baryogenesis between the effective temperatures of 30 GeV and 
0.2 GeV, if the only helds that are allowed to participate in the short distance physics are 
those of the Standard Model. The upper bound on the temperature is motivated from the 
cosmological bound on KK graviton decay while the lower bound is purely for our calcu- 
lational feasibility (related to the QGD phase transition). The time-dependent suppression 
coefficient was motivated from the generic possibility that with the Arkani-Hamed-Schmaltz 
mechanism of suppressing proton decay, the proton decay suppression factor can be time- 
dependent. 

Unfortunately, if the fundamental scale setting the size of the unsuppressed dimension- 
6 operator is forced up to values Mj > 100 TeV by the bound from neutron-antineutron 
oscillations, the branching ratio of 5-violating interactions to 5-conserving interactions is 
set by the quantity {180 rmGeV/Mf)‘^ < 10“^^. We have shown this to be too small to 
allow sufficient baryogenesis for most cases. Another way to view this bound is that as the 
dimension-6 operator creates baryon asymmetry, the same reactants participate in creating 
entropy through 5-conserving operators. The competition between the two types of reac¬ 
tions is almost always dominated by the 5-conserving one. Hence, within the perturbative 
setting where the 4D effective action is valid, sufficient baryogenesis in this scenario is in 
most cases impossible almost independently of cosmology. 

There is one scenario which cannot simply be ruled out by the general analysis. This in¬ 
volves out of equilibrium c quarks interacting with equilibrium b quarks from a temperature 
of about 30 GeV until a temperature of 0.2 GeV. The reason why this scenario evades the sup¬ 
pression is that the primary entropy generation channel is through the b-b self-annihilation 
into c-c instead of the b-c weak interaction, since iWfep ~ 10“^. Furthermore, given that 
nonstandard cosmology of braneworlds may give a slow expansion rate, one can have baryo¬ 
genesis persist for a long time period, giving a large integrated pileup. However, given that 
c quarks are produced efficiently both by the gluons that keep the b quarks in equilibrium as 
well as the b quarks which undergo decay, setting up the necessary hierarchy in densities to 
keep c out of equilibrium while b in equilibrium seems difficult. However, this scenario does 
pass the general tests not relying upon the specihcs of the out of equilibrium mechanism. 

Since neutrinos are almost massless and have vanishingly small tree level annihilation 
cross section in the limit of small y/s, one would naively think that this channel might be 
weak enough for the baryon number-violating operator to compete with it. Unfortunately, 



this turns out not to be the case. 

It is interesting to note that even with just the Standard Model helds, the presence 
of nonrenormalizable operators implies the existence of additional CP-violating rephasing 
invariants beyond the Jarlskog parameter. In any theory, the number of physical CP- 
violating phases is given by the number of couplings in the Lagrangian that are allowed to 
be complex, minus the dimension of the symmetry group which describes phase redefinitions 
of the helds: introducing new effective operators into the SM automatically creates more 
CP-violating invariants, even if all the couplings of the new operators are real in some 
basis. Hence, it is not the smallness of the Jarlskog invariant that prevents the success of 
our baryogenesis scenario. It is merely the fact that unless neutron-antineutron oscillations 
are suppressed by some means other than the fundamental scale, the phenomenologically 
acceptable scale of 100 TeV is still too large for dimension-6 operators to play a signihcant 
role for baryogenesis below the temperature of the electroweak phase transition. 

As with any “no go” arguments, there are many loopholes in our conclusion (besides the 
one remotely possible scenario that we already mentioned). First, we neglected any possible 
nonperturbative physics. The 4D coupling constants are inversely proportional to the volume 
of the extra dimensions in which the gauge helds propagate. Thus, if the volume of the extra 
dimensions was small at some epoch in the early Universe, then nonperturbative physics may 
dominate (see Appendix given that the perturbative contributions to ub/s are insufficient 
in our scenario. For example, the SU{2)l instanton effects may become unsuppressed, giving 
rise to a scenario similar to that of sphaleron transitions at the electroweak phase transition. 
Also, for effective temperatures below about 0.2 GeV, non-perturbative QCD effects will be 
relevant, since the quark degrees of freedom become conhned. Although this is not likely 
to change the suppression of the dimension-6 operator, the kinematics and matrix elements 
will be very different from our perturbative calculation done in the regime of deconhned 
quarks. 

Secondly, we have neglected all effects of the KK modes; or, in the higher dimensional 
picture, we have neglected the fact that the space may be inhomogeneous in the higher 
dimensions. For example, one may envisage a scenario in which there is a phase transition 
which localizes the wave functions of the quarks and the leptons. A proper description of 
quark-lepton separation in higher dimension will require analysis beyond the zero mode. 
However, this type of effect will most likely play a role in enhancing the nonequilibrium 
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condition rather than changing the baryon number-violating branching ratio. Hence, our 
conclusion is most likely robust with respect to this assumption. 

Thirdly, we have neglected all CPT violation effects that must exist because of the time 
dependence of the quark-lepton separation. If the time scale associated with the quark- 
lepton separation process is very short, then there may be significant contributions from 
these effects to enhance the baryon number violating channels. However, such models will 
probably be severely constrained by the restriction on bulk graviton production, just as the 
reheating temperature is severely constrained. 

Hence, although there are loopholes and caveats to our “no go” claims, it seems fair 
to conclude that with only perturbative physics, degrees of freedom beyond the Standard 
Model helds must play a signihcant role in baryogenesis if the effective temperature of 
the Universe never exceeds 30 GeV (which is a very conservative upper bound for the 
maximum temperature in models with large extra dimensions, if the fundamental scale is 
to be accessible to collider experiments). The requirement of beyond the Standard Model 
held content is not an obvious statement, considering that we hardly impose any restrictions 
on the cosmology and that there is an unsuppressed dimension-6 operator as well as an 
effectively large CP violating phase available for baryon asymmetry generation. Of course, 
our “no go” claims would be signihcantly relaxed if the dimension-9 operator responsible 
for neutron-antineutron oscillations is suppressed by some symmetry mechanism which is 
introduced by hand to supplement the geometrical suppression of proton decay operators. 

In closing, note that even if the fundamental scale can be lowered by evading the bound 
from neutron-antineutron oscillation, there is a signihcant challenge in this type of scenario 
as the quarks (and/or leptons) must be forced out of equilibrium at rather high temperatures. 
To accomplish this, it is possible to lower the Planck scale during baryogenesis by having a 
small extra dimensional volume during that period. However, there are severe restrictions 
coming from bulk graviton production in such cases. Furthermore, it is an extremely difficult 
challenge to hnd a time dependent potential for the scalar helds localizing the quarks and 
the leptons, such that the initial quark-lepton separation departs sufficiently from today’s 
equilibrium value in a natural manner. We leave variations on our model that may generate 
sufficient baryon asymmetry to future studies. 
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APPENDIX A: ADEQUATE SUPPRESSION OF PROTON DECAY 

To hnd the minimum value of L consistent with proton stability, we consider the decay 
process in the four-dimensional low energy theory, where we have argued that four-fermion 
operators give the largest contribution. (Processes that cannot be written as four-fermion 
effective operators, for example wormholes, are negligible for the value of L that we hnd is 
needed.) The proton decay rate calculation for these operators is analogous to muon decay 
in the Fermi theory, giving an estimate of 

Tp K, (Al) 

for the “partial lifetime” of decay into any allowed channel, where A is a dimensionless 
constant of order 1. Mode-dependent experimental lower limits on the proton lifetime 0] 
require Tp > 2 x 10^^ years for some modes, which would translate into Mf > 3 x 10^® GeV 
in the absence of the suppression factor Vg- Thus 

a < (M/ ■ 3 X 10^® GeV)2 ^ iq2i-23 _ ^48-53 (^^2) 

where the range quoted is for a fundamental scale ranging from 100 to 10 TeV. The conclusion 
is that about 50 e-foldings of suppression are required at the present epoch. This is easy to 
obtain if Vg ~ e~^ ^ and /i is not too small relative to Mf, but is marginally consistent with 
perturbativity of SU(2)vk, which begins to break down at MjR^ > 30 (recall that > L), 
for proton decay mediated by states of mass Mf, for which a ^ e with k of order 1. 


APPENDIX B: OTHER POTENTIAL SOURCES OF B VIOLATION 


One should also consider nonperturbative processes which may affect B violation and 
baryogenesis [^, Even without the introduction of R-violating operators, the baryon 
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number current is only conserved up to a SU(2)vi/ anomaly, and it is well-known that topo¬ 
logically nontrivial gauge configurations can lead to a change m. B + L, either by quantum 
tunneling via instantons, or thermal excitation via sphalerons 1^^. The rates for B vio¬ 
lation by instantons and sphalerons also have exponential suppression factors and 

respectively, where Ms is the sphaleron energy given by Ms ~ 2M\Y/aw |^- If the 
total length of the fifth dimension were to change with time then, assuming a constant 
five-dimensional gauge coupling these processes could become important. Given the re¬ 
lationship gl = g1 /for the weak gauge coupling, instantons would be exponentially larger 
for smaller R^] the dependence is just g-const.xRs^ instanton enhancement at small 

would shadow the increase in processes mediated by massive states. However, it can 
easily be checked that the instanton rate today is many orders of magnitude smaller than 
the bound on such perturbative processes, therefore instantons can never have dominated 
over the perturbative operators, unless the perturbative operators were further suppressed 
by many orders of magnitude (for example, if they vanished by an exact discrete symmetry). 
The sphaleron rate is more difficult to estimate, since in addition to the temperature depen¬ 
dence one would have to model how the Higgs v.e.v. changes in response to changing R^. On 
the simplest assumption of a constant v.e.v. the sphaleron rate is actually independent of R^ 
and, as in four dimensions, depends just on T. This is a counterexample to the claim that 
all processes are exponentially suppressed with increasing i?; it is understood by noting that 
the formation of sphaleron configurations is a classical, thermally-activated process, in con¬ 
trast to quantum effects which are suppressed by a massive propagator or by the Euclidean 
action. Unsuppressed sphaleron transitions “wash out” any baryon asymmetry generated at 
temperatures above the electroweak scale on the assumption of a very low maximum 
temperature of the Universe, which is required in low-scale models, sphaleron effects will be 
completely negligible. 


APPENDIX C: BOLTZMANN EQUATIONS AND ENTROPY 

Consider first the entropy of one species of particles. Divide the phase space of particles 
into cells and let the cells be numbered i = 1,2,3,... where each cell contains iVj identical 
particles. Let Gi be the phase space volume that can be occupied by one particle in the cell 
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i. The total phase space of the i’th cell is 


AT, = 


G 


Ni 


m 


(Cl) 


if one assumes 1 iVj Gj and negligibly correlated distributions. The total phase space 
volume is 

Ar = nAr, (C2) 

i 

and the dehnition of entropy is 5* = In AT. Using the approximation IniV! ^ iVln(A)^ we 
can express the entropy as 

S = 


d^pd^x 


(27r) 


f{p,x) In 


fiP:X) 


where 


n(x = 


d?p 

{27rf 


f{p,x) 


gives the number density of particles. Hence, one can dehne the entropy current 


= 


d^p 

{2nY 


f{p,x) In 




f{p,x)\ pO 


(C3) 


(C4) 


where we dehne the entropy density as usual as 




(C5) 


Now, using the metric 

= dt^ — a^{t)ddY 

one has 

s^.,. = i + ?>Hs 

and hence the entropy density Boltzmann equation becomes 

f i 


(C6) 

(C7) 


(C8) 


where G/ corresponds to ith collision operator for the particle species / in the Boltzmann 
equation and dT represents the appropriate momentum phase space integral. Neglecting 
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CP violations which are not important for the total entropy generation, we can write the 
collision terms in terms of S-matrix amplitudes as 


-J 2 J 2 j dTC([ff]\n[ff] = dr( 27 r)^( 5 (^) (momenta) x 

/ * i 

(k.h, - /„/*)|M(c., rf.; ai, 6.)P + 

JciJdi 

J (ir( 27 r)^( 5 *'^ymomenta) x 

i 

(fe. - U,H)\M(g., ft.; e.)|= 11117^] + 

9i J 

J (ir( 27 r)"‘( 5 ‘'^ymomenta) x 

i 

/• • • 1 n 1 \i7..ri2i riiiifial states, 

(initial states — hnal states) |iW| ln[—— - J + ... 

hnal states 

+AC (C9) 


where the phase space integral 


rfr^n 


g d^p 


(CIO) 


(27r)3 2E 

is over all particles and momenta, and AC corresponds to the entropy change produced by 
the decay of a classical held conhguration. If a zero-temperature classical held conhguration 
(such as an inhaton vev) decays to particles Xi and X 2 and the Bogoliubov coefficient for 
the particle production is \(3k\‘^ (where the physical momentum is k/a), then we have the 


ehective collision operator 


AC = IdTCf'lh,]\n[h,]+1dTCplfx,]Hh,] = - (Cll) 


contributing to the entropy density Boltzmann equation. 

With our expression for the entropy, one can prove the following little theorem: If baryons 
are generated through particle / which is out of equilibrium undergoing arbitrary reactions, 
then 

ub < —T [ dtT-BY{t)(i^nf, (C12) 

O'* J 


and 


where 



rBv(t) 

. Tt it). 


p = ~ 0 . 02 ^, 

o' ^ 


(C13) 

(C14) 
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T is the effective energy scale which characterizes the kinetic energy of the on shell particles 
in the loop / is the loop integral suppression factor, v is the Higgs VEV of about 180 
GeV, and Fbv is the tree-level baryon number violating reaction rate. The maximum in 
Eq. ( |G10| ) is taken during when most of the baryogenesis is taking place, and Ft corresponds 
to the reaction rate of / which does not violate baryon number. 

The argument validating this claim goes as follows. Starting from Eq. m , one can 
use the momentum conserving ^-function and the form of the equilibrium distribution to 
obtain Eq. (|C12|) if one assumes that the baryon number-violating operator coupling is of 
order 1 (maximal). Now, consider the / evolution. Defining Nf = UfO?, we can write the 
H-conserving part of the Boltzmann equation as 


Nf = -J2 + ACa^ + ... 


(CIS) 


where F* are baryon number-conserving reaction rates, are equilibrium number densities 
per comoving volume, Uj are the ratios of final state number densities to initial state equi¬ 
librium number densities that usually appear in the Boltzmann equations, and ACa^ is the 
external source for /. Since we are assuming maximal out of equilibrium, the equilibrium 
terms involving reverse reactions can be neglected. Hence, we can write 


Nf ^ -Tt{t)Nf + ACa^ 


(C16) 


which governs Nf neglecting the baryon number violating terms. Hence, we can write 

tit') 

Furthermore, we can write 


ns < - Nft')) 


(C17) 


ns < 


2ri 


BV 


dt'AC{t')a^{t') - Nf{t) + Nf{ti) 


■p 

t -I max 

From Eq. (|C5|), we see that 

sa^ > Nf{ti) + f dt'AC{t')a^{t') 

J 

Hence, we find Eq. ( P13|) , as promised. 


(C18) 


(C19) 
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APPENDIX D: SCALING OF TEMPERATURE 


Here we discuss the physical range of z/ in T oc a The generalized hrst law of thermo¬ 
dynamics tells us 

dE = -PdV + TdS + Jpdt (Dl) 

where Jp is the external source/sink of energy for the photons and for other relativistic fluids 
in thermal equilibrium with the photons. Here we have separated the 3 dimensional pressure 
P from the component associated with the extra dimension, whose contribution to dE we 
write as Jpdt. If we assume a constant “effective” average equation of state P/p ^ uj, set 
dS = 0, and assume that dV = a^, then we hnd, trivially, 


P 

— = exp 
pi 


^dt] I - 

a-^p , 


—3(u;+l) 


(D2) 


For photons and other relativistic particles in thermal equilibrium conhned to 3 spatial 
dimensions, the equation of state is a; ~ 1/3 and p oc T^. This implies 


T = Ti exp 


^dt 

a-^p j 


-1 


(D3) 


which implies that v = 1 unless the source term / dt is signihcant. If Jp > 0, then z/ < 1, 
otherwise z/ > 1 . 

An example of a situation in which z/ < 1 is when the bulk helds are in thermal equilibrium 
with the brane helds through 

Jpdt ^ —Pi+d a^dz'^ (D4) 

where P 5 is the bulk pressure and ^ is the scale factor for the extra dimension. 
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